Staggered- Vortex Superfluid of Ultracold Bosons in an Optical Lattice 
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We show that the dynamics of cold bosonic atoms in a two-dimensional square optical lattice pro- 
duced by a bichromatic light-shift potential is described by a Bose-Hubbard model with an additional 
effective staggered magnetic field. In addition to the known uniform superfluid and Mott insulating 
phases, the zero-temperature phase diagram exhibits a novel kind of finite-momentum superfluid 
phase, characterized by a quantized staggered rotational flux. An extension for fermionic atoms 
leads to an anisotropic Dirac spectrum, which is relevant to graphene and high-T c superconductors. 



The experimental realization of ultracold atomic gases 
loaded into optical lattices has opened up a unique path- 
way to study quantum phase transitions of many-body 
systems [lj . Exploiting the rich internal atomic structure 
and the great versatility in the engineering of optical po- 
tentials, increasingly complex optical lattice models have 
been proposed for which unconventional quantum 

phases are predicted. Recently, it has been suggested 
that optical lattices can be used to simulate the quan- 
tum behavior of charged particles in a two-dimensional 
(2D) lattice subjected to a homogeneous magnetic field 
[J]. This system is known to exhibit a wealth of inter- 
esting physics, such as the famous Hofstadter butterfly 
single-particle spectrum Q or the integer and fractional 
quantum Hall effects [f|. Much less is known, however, 
about charged particles moving in a 2D lattice subjected 
to a staggered magnetic field. In a pioneering work, Hal- 
dane has shown that the integer quantum Hall effect may 
occur in this system as a result of broken time-reversal 
symmetry Q. Concerning the single-particle spectrum, 
only recently numerical studies have revealed the con- 
nection to the Hofstadter butterfly [§]. The technical 
difficulty to engineer magnetic fields alternating on the 
spatial scale of condensed matter lattices has constrained 
experimental studies to magnetic fields modulated on a 
mesoscopic scale [§]. 

In this letter we show that a staggered magnetic field 
in a lattice can be realized for ultracold bosonic and 
fermionic atoms in a 2D optical lattice. We then present 
a mean field theory for the bosonic system and construct 
the zero-temperature phase diagram shown in Fig. [T] 
For small magnetic fields, the spatially uniform (k — 0) 
superfluid phase and the Mott insulating state, known 
from the conventional Bose-Hubbard model, are repro- 
duced. When the magnetic flux per plaquette <E> exceeds 
$o/2 (see Fig. a) for the definition of plaquette), where 
$0 = hc/e is the fundamental flux quantum, we find that 
a novel finite momentum (k = n) superfluid phase is real- 
ized. This superfluid phase is spatially modulated and is 
characterized by quantized fluxes of alternating sign for 
adjacent plaquettes. Our work thus points out the pos- 
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FIG. 1: (color online). Phase diagram with respect to the 
chemical potential fj,, the interaction parameter U, the hop- 
ping amplitude J, and the scaled magnetic flux W. Within 
the three-dimensional lobes a gapped Mott insulator phase 
prevails. Outside these lobes the system is superfluid. For 
small magnetic fields (W/J < 1) the superfluid is spatially 
uniform, while for large magnetic fields (W/J > 1), a quan- 
tized staggered rotational flux arises. The white dashed line 
indicates the (fJ,/U = 2 — \/2) plane analyzed in detail in 
Fig.0 



sibility to realize a novel spatially modulated superfluid 
with ultracold bosons, which exhibits a characteristic mo- 
mentum spectrum. An extension of our work for fermions 
offers the possibility to simulate various strongly corre- 
lated systems, such as the mean-field Hamiltonian of Af- 
fleck and Marsto n llQfl , proposed in the context of high-T c 
superconductors [llj, and the case of massless Dirac spec- 
tra, as realized in graphene [lij ]. 

Effective Hamiltonian.~We first consider bosonic 
atoms trapped in a 2D square optical lattice potential 
Vl(y) = — Vo(s\n 2 (kx) + sin 2 (ky)) with an additional 
time-dependent potential Vr(t, t) = kVl(f) cos(2S'(r) — 
ilt), S(r) = tan _1 [(sin(fca;) — sin(fcy))/(sin(fca;) + 
sin(fcy))], which acts as a collection of micro-rotors, each 
applying angular momentum to a single plaquette, with 
alternating signs for adjacent plaquettes. Here, k = 2tt/X 
and A is the wavelength of the optical potential. The term 
Vft(r, t) can be implemented experimentally by means of 
a bichromatic light-shift potential fl3| . The well depth 
Vq > 0, the oscillation frequency f2 and the coupling 
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FIG. 2: (color online), (a) The bipartite lattice comprises 
sites labelled by "A" and "B." The gray rectangle identifies 
a single plaquette. The solid arrows indicate the tunnelling 
currents driven by Vr(t, i). (b) Phase diagram within the 
(n/U = 2 — \/2)-plane spanned by the dashed white line and 
the W axis in Fig. [T] 



strength k G [0, 1] arc adjustable parameters in experi- 
ments. We denote the sublattices of the bipartite square 
lattice by A and B (cf. Fig. [2^a)) and define four vec- 
tors ei = — e 3 = (A/2)x, e 2 = — e 4 = (X/2)y connecting 
each A site to its four B neighbors. Upon the assump- 
tion that the atoms in the optical lattice are restricted to 
the lowest Bloch band and that Vr(t, t) does not induce 
interband transitions, the system can be described by a 
Bose-Hubbard model with time- varying hopping and en- 
ergy offset [3, HI 



H{t) = - J2 K + H - c -} 

rSA,i=l-4 

+ e r (t) h T + 2 U J2 n r (n r - 1) , (1) 



r<£A<SB 



where a r and a£ are the boson annihilation and cre- 
ation operators on site r obeying the canonical com- 
mutation relation [a r> ttj/] = o~ r y, n T = aj a r is the 
number operator, Jj(t) = J + (— 1)'kVoXi sin(Of) de- 
notes the anisotropic time- varying hopping, e r eA,B(t) = 
±2kVqX2 cos(£2t) is a time-varying energy offset, and 
U is the onsite interaction energy. In terms of the 
Wannier function of the lowest band w(r), we have 
J = J dxdy w* (x + A/4, y) [- £ V 2 + V L (r)]w(x - A/4, y) , 
Xi = Jdxdyw*(x + A/4, y) [sin 2 (kx) — cos 2 (ky)]w(x — 
A/4, y), X2 — J dxdy \w(x,y)\ 2 [cos(kx) cos(ky)] and U oc 
(a s /m) J dxdy \w(x, y)| 4 , with the atomic mass m and 
the s-wave scattering length a s . We note that by us- 
ing Feshbach resonances, the s-wave interaction can be 
considered as another experimentally tunable parame- 
ter. Upon expressing all terms with harmonic time- 
dependence in terms of a quantized auxiliary bosonic field 
and integrating out this field, similarly as in Ref. 0], we 
find the effective Hamiltonian 



H, 



eg ' 



J2 {|c|e i0 (- 1 ) 1 ata r+e! +H.c.} 

+ hj n r (n r - 1), (2) 



reA,Z=l-4 



where the anisotropic complex hopping amplitude is 
given by |c| = V J 2 + W 2 , 9 = taxT x {W/J) with W = 
2k 2 V^xiX2/^- The Aharonov-Bohm phase 6, picked up 
by the bosons when tunnelling along the edge of a pla- 
quette, may be interpreted as resulting from an effective 
staggered magnetic field with a magnetic flux per plaque- 
tte $ = (2f3/7r)<E>o- This artificial gauge field arises due 
to the 7r/2 phase lag between the time- varying hopping 
terms Ji{t) and the energy offset terms e rE x,B(i) ; which 
assigns a sense of rotation to each plaquette. In Eq. ([2|) 
we have neglected non-local negative ring exchange terms 
with energy (kVqXi) 2 /2hil which may be tuned to be 2 
orders of magnitude smaller than all other terms. We also 
omitted a nonlocal term composed of summands n r n r ', 
describing interaction between distant lattice sites medi- 
ated by the retroaction of the atoms upon the light-shift 
potential. This term is expected to be significant, if the 
light-shift potential is produced inside a high-finesse op- 
tical cavity 14j. For conventional light-shift potentials 



produced by superposition of laser beams, which we con- 
sider here, it is irrelevant. Nonetheless, this term does 
not alter the mean-field results presented below. 

For an optical lattice loaded with fermionic atoms, a 
similar derivation leads to an effective staggered mag- 
netic field as well. The single-particle spectrum can 
be obtained for both, fermions and bosons, by ex- 
pressing the kinetic terms in Hamiltonian ^ in mo- 
mentum space with two independent amplitudes ak = 
SrgA a r e zk r , 6k = E re s a r e lk r , and applying the canon- 
ical transformations dk = [(e k /l e k|)ak + bk]/y/2, /3k = 
[-(Sk/MR + k]/V2. Here, e k = 2\c\e iB cos[fc+a] + 
2\c\e~ lS cos[fc~a] is the generalized lattice dispersion, 
with fc ± = (k x ± ky)/2 and a = X/y/2. We then find 

E± = ±2|c|[cos 2 (fc+a) +cos 2 (fc"a) 

+2 cos (fc+a) cos (AT a) cos (26*)] 1 / 2 , (3) 

where the lower (upper) branch corresponds to the in- 
(out-of-)phase mode /3k (dk) between the sublattices. For 
the fermionic system at half-filling, the single-particle 
spectrum exhibits two incquivalcnt anisotropic Dirac 
cones; see Fig. [3] Their slope can be tuned via the 
staggered magnetic field, thus resembling the physics of 
graphene with anisotropic hopping arising under uniax- 
ial pressure 15(. Furthermore, at 8 = ir/4 and negligible 
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interactions, the system simulates the Affieck-Marston 
Hamiltonian [Toj . wherein a staggered 7r-flux phase is 
proposed to describe the pseudogap regime of high-T c su- 
perconductors. Although several experiments were sug- 
gested to probe the relevance of this phase, the problem 
remained open due to unavoidable effects of disorder [lj| . 
The system we consider thus offers an opportunity to 
study such phases in a highly controllable environment. 

Mean-field theory for the bosonic system.-We antici- 
pate the bosonic system to exhibit the well-known quan- 
tum phase transition from a superfluid to a Mott insulat- 
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|*o> « (ft) N \0) = (EreA®B&l) N \Q)- For 8 > tt/4, new 
absolute minima develop at the Brillouin zone edges k = 
(±7r/a, ±7r/a) = tt for which condensation takes place; 
see Fig. [3] Because of the equivalence of the four minima 
in the reciprocal space, the new many-body ground state 
can be written as cx ) N \0) = (J2 r (°r + i ^t+e 2 - 
&t +ei+e2 - iat +e J) N \Q) with r = 2r' and r' £ A © B. 
The angular phases of the order parameter differ by 7r/2 
between neighboring lattice points, and there is a quan- 
tized flux on each plaquette, with alternating sign for ad- 
jacent plaquettes. Thus, for the BEC at k = n, the sys- 
tem is characterized by a vortex-antivortex lattice with a 
periodicity X/V2, namely, a staggered- vortex superfluid 
phase. This phase possesses a definite chirality that is 
commensurate with the staggered magnetic field, a fea- 
ture which is not present in the uniform superfluid phase. 
To confirm that distinct superfluids are stabilized, we em- 
ploy a variational mean-field ansatz for the ground state 
\{ } o) = (e" l «/ 2 cob(<t)PI + e*«/ 2 sin(CT)/3t)^|0) and mini- 
mize the expectation value with respect to the Hamilto- 
nian @. The uniform superfluid (a = <jq = 0) and the 
staggered- vortex phase (a = ao = tt/2) are indeed the 
absolute minima of the mean-field energy for 8 < 7r/4 
and 8 > tt/4, respectively. Furthermore, the stability 
of both phases can be verified by examining the en- 
ergy cost for deviating from the condensate (H)\ ao ± s \ = 
Emf + eHN\c\ [V2| sin(7r/4 - 8)\ + (UN)/(8\c\)] +0(e i ) 
with cT-o € {0, tt/2} and Emf = -4N\c\ cos(0). Finally, 
we note that as the system is tuned across the 8 = 7r/4 
line, ao changes discontinuously by a value of 7r/2, sug- 
gesting that the two superfluid phases are separated by 
a quantum first-order phase transition line within this 
variational mean-field analysis. At this line, both phases 
are degenerate in their mean-field energies, with a finite 
energy barrier A ~ UN 2 /8N S between the two minima. 

Fluctuations. -Following the Bogoliubov theory for a 
weakly interacting Bose gas, we derive the energy spectra 
of the superfluid phases. We consider the Hamiltonian 
@ in the grand-canonical ensemble and make the sub- 
stitution for the condensation mode /3 ko — > + Pk > 
where Nq is the condensate number and ko = (ko = 
7r) for the uniform superfluid (staggered-vortex) phase. 
Choosing the chemical potential at its mean-field value, 
fi = — | ejc n | + noli/2 and keeping terms up to quadratic 
order in the fluctuations, the Hamiltonian becomes, 




FIG. 3: (color online). Single-particle spectrum for (a) 6 < 
7r/4, with the minimum at the origin, and (b) 9 > tt/4, with 
minima at the Brillouin zone edges (±7r/a, ±7r/o). 

ing phase as the parameter C//4|c| is increased The 
phase boundary between the two phases is determined by 
a straightforward generalization of the functional integra- 
tion method presented in Ref. 17 1. The partition func- 



tion is written in terms of an imaginary-time functional 
integral Z = jT>a*T>aexp{—S[a* ,a)/h} with the action 

S[a*,a] = / S/fcflT * E r <(r) (M T - ^(r) + 
where T is the temperature and /i the chemical poten- 
tial. In the Mott regime, we introduce the order parame- 
ter ip r (t) (a Hubbard-Stratonovich field) to decouple the 
hopping term and upon integrating out the boson fields 
(a*, a), the effective action (up to quadratic order) gives 
the zero-temperature quasiparticlc(holc) energy disper- 



sion, 



qp,qh 



-M+f(2n-l)-^±^ k , (4) 



where huj^ = ^/ 1 e k | 2 — (An + 2)|e k |C/ + U 2 is the energy 
required for creating a quasiparticle-quasihole pair and 
n is the lattice filling factor, given by the ratio between 
the number of bosons N and the number of sites A^ s . 
The phase boundary between the Mott and the super- 
fluid states can then be determined by the condition that 
the excitation becomes gapless tuo-^ — 0, which is shown 
in Fig. [T] as the generalized Mott lobes. For brevity we 
henceforth concentrate on the (fi/U = 2 — v / 2)-plane 
spanned by the dashed white line and the W-axis in Fig.[T] 
[cf. Fig. GJb)]. We note that for 6 > vr/4 the low en- 
ergy excitations in the gapped Mott phase carry a finite, 
rather than zero, lattice momentum. This implies that 
the associated critical phenomena at the phase bound- 
ary between the Mott-insulator and the staggered vortex 
phase could be different from the usual universality class 
of the 0(2) quantum rotor model for the superfluid-Mott 
transition in the Bose-Hubbard model [l8j |. 

In the superfluid regime, however, the staggered mag- 
netic field drives the system into distinct superfluid 
phases. Even though interactions induce a finite quan- 
tum depletion, for a system of N weakly interacting 
bosons, Bose-Einstein condensation (BEC) takes place 
at the lowest single-particle state. For 8 < tt/4, BEC 
occurs at the k = (0, 0) = state, giving rise to a 
uniform superfluid with the many-body ground state 



-|e k |-|e ko | + -n [/]/3 k /3 k 



in LT^ k>ko (4dl k + /3 k /JL k ) 



-~no^.Bk ) k a k /3l k + H.c. 
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FIG. 4: (color online). Schematic of the momentum spectra 
for the uniform superfluid (a) and the staggered-vortex su- 
perfluid (b). The colored areas in the centres illustrate the 
first Brillouin zones with (b) and without (a) the staggered 
magnetic field. 



where no = Nq/N s is the condensate density, 
^4k,o = -S kj7r = 1 + exp(-2i<£ k ), 5 k ,o = A k<7T 
1 — exp(— 2iipk), and ip^ = arg(e k )- The Hamil- 
tonian can be readily diagonalized by a Bo- 
goliubov transformation to yield the spectrum 



|Ek| + l e k | H 

where G kj o = cos 2 (y k ) 
and Gk,Tr = sin 2 (<pk) - 
The lower branch of 



n C/|e ko | ± n U\ek\y/Gk,k a , 

+ 4|e |[|e |/noC/ + l}/n U 
- A\e„\[\e«\/n Q U + l]/n Q U. 
the spectrum is lin- 



ear and gapless at long wavelength fiwk,k ~ 
yj\c\ cos(6> - fc a/2)[4|c| cos(6» - k a/2) + n U] (k - k ) 
corresponding to the Goldstone mode of the broken 
gauge symmetry. 

Experimental signatures. -The characteristic momen- 
tum spectrum of the staggered-vortex phase provides a 
clear signature to identify this state experimentally by 
imaging momentum space using standard ballistic ex- 
pansion techniques. The momentum distribution can 
be expressed as (tft(k)¥(k)) = \W(k)\ 2 S B (k)S P (k), 
where W(k) is the Fourier transform of the Wan- 
nier function, 5b (k) = I Ereaqb g*2k-R|2 j g ^ e s ^ ruc _ 
ture factor of the Bravais lattice, and 5p(k) = 



Ate{l,2,3,4} 



3 ik- (r 



is the structure factor 



of a plaquette. Here, r„ denote the four corners of a pla- 
quette and a Tv are the corresponding boson operators. 



2 with 
= vir/2 



We may write S P (k) = 1 X^^e{i,2,3,4} e ' 
ip u = for the uniform superfluid \^o) and ipv 
for the staggered vortex superfluid I'Fjr)- As illustrated in 
Fig. [3J the two cases display distinct structures of Bragg 
maxima, directly observable in experiments. 

In conclusion, we have shown that anisotropic and 
time-varying hopping terms in a 2D optical lattice give 
rise to an effective staggered magnetic field. For the 
bosonic system, it leads to a novel kind of superfluid 
phase characterized by a quantized staggered rotational 
flux. For the system realized with fermionic atoms, it 



gives rise to anisotropic Dirac spectra at half-filling. The 
tunability of the interaction terms and the addition of op- 
tical disorder potentials allow for systematic simulations 
of various strongly correlated systems, such as graphene 
and high-T c superconductors. Another exciting direc- 
tion for future work could be the search for quantum 
Hall physics in the system. Finally, we remark that the 
experimentally accessible Hamiltonian offers a wider 
parameter space than presently considered in this work. 
The inclusion of non-negligible ring exchange interac- 
tions, for example, may offer the opportunity to realize 
exotic quantum insulators Q. 
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